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Abstract 

Evolution of a universe with homogeneous extra dimensions is studied with 
the benefit of a well-chosen parameter space that provides a systematic, useful, 
and convenient way for analysis. In this model we find a natural evolution 
pattern that entails not only stable extra dimensions in the radiation-dominated 
era, thereby preserving essential predictions in the standard cosmology, but also 
the present accelerating expansion while satisfying the limit on the variation of 
Newtonian gravitational constant. In this natural evolution pattern the extra 
dimensions tend to be stabilized automatically without resorting to artificial 
mechanisms in both the radiation-dominated and the matter-dominated era, as 
a wonderful feature for building models with extra dimensions. In addition, the 
naturalness of this evolution pattern that guarantees the late-time accelerating 
expansion of a matter-dominated universe presents a solution to the coincidence 
problem: why the accelerating phase starts at the present epoch. The feasibility 
of this evolution pattern for describing our universe is discussed. 
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1 Introduction 



The observations by Supernova Cosmology Project and Supernova Search Team have 
suggested in 1998 that the expansion of the present universe is accelerating [TJ |2]. 
This conclusion is reinforced recently in 2003 by WMAP measurements [3]. One gen- 
eral conclusion from these measurements and the CMB observations in recent years 
in El 0] is that the universe has the critical density, consisting of 1/3 of ordinary 
matter and 2/3 of dark energy with a negative pressure [H] (such that Px/Px < —0.78 
[3]). Although this acceleration may be driven by the existence of a positive cos- 
mological constant (vacuum energy) [7], there remain other interpretations of the 
accelerating expansion, such as "quintessence" (a slowly evolving scalar field) [HI E] 
and the presence of extra dimensions [TUl ITT] . 

The presence of extra dimensions is required in various theories beyond the stan- 
dard model of particle physics, especially in the theories for unifying gravity and other 
forces, such as superstring theory. Extra dimensions should be "hidden" (or "dark" ) 
for consistency with observations. Various scenarios for "hidden" extra dimensions 
have been proposed, for example, a brane world with large compact extra dimen- 
sions in factorizable geometry proposed by Arkani-Hamed et al. ^2] (see also jl3j). 
and a brane world with noncompact extra dimensions in nonfactorizable geometry 
proposed by Randall and Sundrum In this paper we will employ the simplest 
scenario: small compact extra dimensions in factorizable geometry, as introduced in 
the Kaluza-Klein (KK) theories [T^ . 

The possibility of generating the accelerating expansion of the present matter- 
dominated universe via the evolution of homogeneous and isotropic extra dimensions 
is first pointed out by Gu and Hwang [TU]. Furthermore, the general idea of unifying 
dark energy sources (i.e. dark matter and dark energy) and dark geometry (e.g. extra 
dimensions) into one has been sketched by Gu fT], who pointed out that dark geome- 
try (instead of dark energy) can be an intriguing candidate for generating accelerating 
expansion. As indicated by Einstein's general relativity, mass (energy) and geometry 
are two faces of one single nature, therefore it is biased to consider only mass (energy) 
has dark part, while believing blindly that all geometry is totally "visible" to our poor 
eyes [TT]. 

In this simple scenario making use of a highly symmetric extra space, there are 
much fewer free parameters (only two additional degrees of freedom, the expansion 
rate and the curvature of the extra space, in addition to those in the standard cos- 
mological models without dark energy) so that it is much easier to be ruled out, 
compared with the quintessence models, by constraints from observations. In par- 
ticular, an essential difficulty of this model ^H] stems from the constraint on the 
variation of the Newtonian gravitational constant, which will be produced along with 
the evolution of extra dimensions that is the key element for generating accelerating 
expansion. 

In this paper we will study a more general case in which the evolution of our uni- 
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verse in various eras, especially the present accelerating expansion era, is governed not 
only by the matter contents (excluding dark energy) therein, but also by the curvature 
of the ordinary 3-space and the evolution of extra dimensions. We will explore the 
feasibility of this model for generating accelerating expansion while satisfying obser- 
vational constraints, especially the variation of the Newtonian gravitational constant. 
Through these studies we will also present several nice features of this model, such 
as the automatic stabilization of extra dimensions and the solution to the cosmic 
coincidence problem — why the energy densities of dark energy and dark matter are 
comparable now (i.e. "why now" problem), or, more precisely (if dark energy is not 
a necessary ingredient for the accelerating expansion), why the accelerating phase of 
our universe starts at the present epoch. 



2 The Higher-dimensional World: Basics 

We consider a {3 + n + l)-dimensional space-time where n is the number of extra 
spatial dimensions. Assuming that both the three-dimensional ordinary space and 
the n-dimensional extra space are homogeneous and isotropic, we use two spatial 
parts of the Robertson- Walker metric to describe this space-time as follows: 

ds^ = de - a\t) + rldnl)^ - h\t) + rldni) , (1) 

where a{t) and h{t) are scale factors, and ka and ki, relate to curvatures of the ordinary 
3-space and the extra space, respectively. The value of is set to be within the 
interval [0, 1) corresponding to the finite size of extra dimensions. Assuming that the 
matter content in this higher- dimensional space is a perfect fluid with the energy- 
momentum tensor 

T'^f3 = diag{p,-pa,...,-Pb,...), (2) 
we can write the Einstein equations (TUl E] as 
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where p is and the energy density in the higher-dimensional world, and pa and pb are 
the pressures in the ordinary 3-space and the extra space, respectively. 

The terms involving the scale factor h in Eqs. Q and (0} are additional terms 
and Eq. (jSJ is an additional equation coming from extra dimensions. As pointed out 
by Gu in if extra dimensions exist but we are so biased that we admit only 
mass/energy can have dark part while blindly believing there is no dark geometry 
(such as extra dimensions) at all, these additional terms from extra dimensions will 
automatically be moved from the left-hand side of the Einstein equations (describing 
geometry) to the right-hand side (describing energy contents) and be treated as some 
sort of effective dark energy sources, as indicated in the following equations [from 
rearranging Eqs. Q^©]: 



Note that an equation such as Eq. (0) or Eq. (jH)) is absent in the standard cosmology, 
that is, a no-dark-geometry believer living in a (3 + l)-dimensional universe will 
never notice this equation, while a higher-dimensional cosmologist will realize it as 
a constraint equation for the behavior of the scale factor h and the pressure in the 
extra space ph- 

Using equations of state pa = Wap and pb = Wbp, we rearrange Eqs. (H)) and (0) 
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(12) 
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Before performing detailed numerical analyses, we can extract essential features 
of the above equations and possible evolution patterns governed by them via new 
variables 



u(t) = a/ a. 
v{t) = 6/6, 

and the parameter space expanded by 
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The use of the parameter space for analyzing the evolution of the universe with extra 
dimensions was first suggested by Gu and Hwang in [TIH lllj. In terms of the above 
new variables, Eqs. Q, (fT^. and (fTSj) can be rewritten as follows. 
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The above three equations indicate that the essential quantities Gp, a/ a and h/h, in 
unit of V?, are all functions of Xa-, Xb and Y for given values of n, Wa and Wb- That 
is, Xai Xb and Y form a set of variables that can specify the state of a universe with 
extra dimensions, provided that the unit of time is set to be the Hubble time at 
each moment. With the use of the parameter space {Xa , Xf, , F) we can plot regions 
for p > 0, p = 0, or p < 0, and regions for a > 0, a = 0, or a < 0, in order to catch a 
basic picture of the state of the universe at each position in the parameter space. In 
addition, drawing flow vectors in the parameter space is useful for sketching possible 
evolution patterns. Flow vectors in the parameter space are vector fields: 
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which, in unit of u, are also functions of {Xa , X;, , F) for given values of {n ,Wai Wb), as 
guaranteed by Eqs. (pi?j) - (PT|) . Accordingly the use of the parameter space {Xa ,Xb,Y) 
makes it possible to study the evolution of the universe before doing brute-force nu- 
merical calculations. 

In the rest of this paper we will only consider the universe with fiat extra di- 
mensions, i.e., kb = 0, for simplicity. Accordingly, in the next section we will make 
use of a two-dimensional parameter space expanded by Xa and Y (i.e. the plane of 
Xb = 0) to study the evolution of the universe in various possible eras, from the 
very early time to the present. [Hereafter we will use X instead of Xa for conve- 
nience.] We will discuss three possible eras: the blazing era, the radiation-dominated 
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era, and the matter-dominated era, following the time order. The existence of the 
blazing era is speculative. In this era the temperature of the universe is so high that 
particles move relativistically in both ordinary dimensions and extra dimensions and 
hence Wa = Wb = 1/ {3 + n) . In the expanding and cooling process of the universe, 
after some time the temperature will be lower than the inverse of the size of extra 
dimensions so that the motion in extra dimensions, i.e. KK modes, cannot be excited 
and particles are basically at the ground state and therefore pressureless in extra 
dimensions thereafter. After that time if the temperature is high enough such that 
dominant particles are still relativistic in ordinary dimensions, we have Wa = 1/3 and 
= 0, as to be designated as "radiation-dominated era" in this paper. Conversely, if 
the temperature is so low in the later time that the universe is dominated by particles 
that are nonrelativistic in both ordinary dimensions and extra dimensions, we have 
Wa = Wb = 0, as to be designated as "matter-dominated era". 



3 The Evolution in Three Eras 
3.1 The Blazing Era 



n=3 v/u u=d/a 




Figure 1: The {v/u ,ka/a^u'^) parameter space describes various states and possible 
evolution patterns of a blazing universe, where the number of extra dimensions is 
specified to be three for demonstration. Flow vectors and four evolution trajectories 
(white dashed curves) are plotted. 

In the blazing era, Wa = Wb = l/{3 + n). In Fig. [T]we plot regions for p > 0, 
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p = or p < 0, which are obtained from Eq. p9p . regions for a>0,a = Oora<0 
from Eq. (j^O)), flow vectors d{X , Y)/{udt) from Eqs. (j2ni)-(j22I) (assuming u > 0), and 
four trajectories from numerical solutions, where the number of extra dimensions n 
has been specified to be three for demonstration. The black area denoted by "p < 0" 
is a forbidden region if the energy density is required to be positive. White dots 
stand for fixed points where flow vectors are zero both in X and Y directions. The 
four evolution trajectories correspond to numerical solutions of Eqs. (jl2p and (jl3j) 
with the initial conditions (kja'^u'^^v/u) = (a) (-0.0625,6.25), (b) (0.0625,0.75), 
(c) (1,-0.9), and (d) (-0.01,-0.1) respectively, representing four possible evolution 
patterns:^ 

(a) acceleration deceleration acceleration, eventually approaching asymp- 
totically the attractor at (—1, 0) with stable extra dimensions and zero acceler- 
ation, possessing negative spatial curvature. 

(b) acceleration deceleration, possessing increasing positive curvature con- 
tribution. 

(c) deceleration — > acceleration, eventually approaching asymptotically the at- 
tractor at |0)~ 3 + y^3{n + 2)/n /(n — 1))| with collapsing extra dimen- 
sions, possessing decreasing positive curvature contribution in the late time. 

(d) deceleration acceleration, eventually approaching asymptotically the at- 
tractor at (—1, 0) with stable extra dimensions and zero acceleration, possessing 
negative spatial curvature. 

It is essential to ask whether in this model the inflation can take place in this 
era without introducing the inflaton field. From Fig. ^ we conclude that the universe 
cannot have a sufficient inflating expansion, which makes the universe nearly flat, 
in the blazing era unless the extra space is collapsing violently or the expansion 
rate of the extra space is much larger than that of the ordinary space (i.e. v ^ u) 
initially. On the contrary, if the inflation occurs before the blazing era such that the 
resultant curvature contribution is extremely small, the universe may evolve nearly 
along the w/w-axis during the blazing era. In this case, for a signiflcant part of initial 



states in the parameter space (more precisely, for v > — 3 + A/3(n + 2)/n /{n — 1) 

if ka/a^u^ ~ as required by the inflation) the universe tends to evolve to the state 
at {kal a^v? ,v /u) = (0,1), that is, the ordinary space and the extra space tend to 
synchronize their expansion rates. This feature provides a natural initial condition, 
ka/a^u^ ~ and v/u ~ 1, for the succeeding radiation-dominated era. 

^Notc that for all the three eras we will only present some among all the possible evolution 
patterns because many others are apparently not suitable for describing our universe. 
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3.2 The Radiation-dominated Era 



In the radiation-dominated era, Wa = 1/3 and Wb = 0. Fig. El shows regions cor- 
responding to different signs of p and a, flow vectors, four fixed points (denoted by 
four white dots), and three evolution trajectories in the parameter space (where n 
is specified to be three for demonstration). The three trajectories are with respect 
to initial conditions (fca/aV,w/u) = (a) (-0.0625,5), (b) (0.0625,5), (c) (1,-0.9), 
representing three possible evolution patterns: 



^ V < V V \ \ ^ A 

V V < v \ \ 

V V V V V \ \ 

V V V V 

V V V v 

»^ V V V 




Figure 2: The {v /u ,ka/ a^v?) parameter space describes various states and possible 
evolution patterns of a radiation-dominated universe, where the number of extra di- 
mensions is specified to be three for demonstration. Flow vectors and three evolution 
trajectories (white dashed curves) are plotted. 



(a) acceleration deceleration, eventually approaching asymptotically the at- 
tractor at (—1, 0) with stable extra dimensions and zero acceleration, possessing 
negative spatial curvature. 

(b) acceleration deceleration, eventually approaching asymptotically the ka/a^v^- 
axis with stable extra dimensions, possessing increasing positive curvature con- 
tribution. 

(c) deceleration acceleration, eventually approaching asymptotically the at- 
tractor at |o,- ?, + ^/^^{nTi)J^ /(^-l))} 

with collapsing extra dimen- 
sions, possessing decreasing positive curvature contribution in the late time. 
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The answer to the question regarding whether the inflation can be generated in 
this radiation-dominated era without introducing the inflaton field is similar to that 
in the blazing era. We conclude from Fig. |21 that the universe cannot have a sufficient 
infiation, which makes the universe fiat enough, in this era unless the extra space is 
collapsing violently or the expansion rate of the extra space is much larger than that 
of the ordinary space (i.e. v ^ u) initially. 

An important feature in this era is that for about half of initial conditions in the 
parameter space the corresponding evolution will eventually approach to the states 
with stable extra dimensions, i.e., v/u — >■ 0. This implies that the extra dimensions 
have a significant possibility to be stabilized automatically without resorting to addi- 
tional artificial mechanisms. For a negligibly small expansion rate of extra dimensions 

we can recover the standard cosmology as follows: With negligible contributions 
from f , zero curvature in the extra space {ki, = 0), and equations of state for this era 
{pa = l/3p and pb = 0), Eqs. (0)^® become 



where Eqs. ()23p and ()24|) are exactly the Einstein equations in the standard cosmology 
and Eq. (^3)) is a redundant equation that can be derived from the other two equations. 
[Indeed the fact that one of the above equations is redundant indicates the existence 
of the solution with v = 0.] Consequently, with small enough v in this model we 
can preserve the essential predictions in the standard cosmology, in particular the 
primordial abundance of light elements produced during Big Bang Nucleosynthesis 
(BBN), in the radiation-dominated era. In addition, this feature provides a natural 
initial condition, ka/a^v? ~ (guaranteed by infiation presumably) and f/w ^ 0, for 
the succeeding matter-dominated era. 



The evolution in the matter-dominated era, especially the possible accelerating ex- 
pansion of the present universe, is the main issue that we will use more length to 
discuss in this paper. In this era the universe is dominated by pressureless matter 
(without introducing dark energy), therefore p = pM and Wa = Wb = 0. For demon- 
stration the number of extra dimensions n is chosen to be three for the plots in the 
parameter space. Fig. El shows regions corresponding to various signs of and a 
and the observational bound to the energy density of the matter (mainly the cold 
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3.3 The Matter-dominated Era 
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v/u 




u=d/a 
v=b/b 
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Figure 3: Various conditions in the {v/u ,ka/a^u^) parameter space for a matter- 
dominated universe are demonstrated. 



u=d/ a 




Figure 4: The {v/u , ka/a'^u'^) parameter space describes various states and possible 
evolution patterns of a matter-dominated universe. Flow vectors and four evolution 
trajectories (white dashed curves) are plotted. 
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dark matter), which is represented by a sohd curve denoted by fl(M) = 0.35.^ Fig. HI 
shows flow vectors, four fixed points (denoted by four white dots), and four evolu- 
tion trajectories in the parameter space. The four trajectories are with respect to 
initial conditions {ka/a^u^,v/u) = (a) (-0.0625,5), (b) (0.0625,5), (c) (1,-0.9), (d) 
(—0.0625, —0.1575), representing four possible evolution patterns (that are somewhat 
similar to those in the blazing era): 

(a) acceleration deceleration acceleration, eventually approaching asymp- 
totically the attractor at (—1, 0) with stable extra dimensions and zero acceler- 
ation, possessing negative spatial curvature. 

(b) acceleration — > deceleration, possessing increasing positive curvature con- 
tribution. 

(c) deceleration acceleration, eventually approaching asymptotically the at- 
tractor at |0)~ 3 + ^/3{n + 2)/n /(n — 1))| with collapsing extra dimen- 
sions, possessing decreasing positive curvature contribution in the late time. 

(d) deceleration — > acceleration, eventually approaching asymptotically the at- 
tractor at (—1, 0) with stable extra dimensions and zero acceleration, possessing 
negative spatial curvature. 

The possible state and the allowed evolution pattern of the universe are con- 
strained by a variety of observations, such as anisotropies of the cosmic microwave 
background (CMB), large-scale structure (LLS), distance measurements of type la 
supernovae (SNe la) etc. In addition to the limit to the variation of the Newtonian 
gravitational constant that entails the a general difficulty for this scenario |16j . 
we will consider constraints on the "acceleration parameter" go = ao/aoMg (i-e. —qo) 
and other essential cosmological parameters such as the pressureless matter density 
parameter Qm and the Hubble constant Hq (i.e. uq) (where the subscript '0' denotes 
the present time). 

As summarized in the bound to \Gn/Gn\ ranges from 10~^^ yr~^ to 10~^° 
yj,-i Jqj^ different observations. In the model with extra dimensions considered here, 
the Newtonian gravitational constant varies along with the variation of the size of 
extra dimensions in the way that 

GN = G/VnOcb-^, (26) 

where Vn is the volume of the extra space. Accordingly, 

Gn/Gn = -nv . (27) 



^The observational bound to the energy density of the pressureless matter is 0.15 ^ fi^/ ^ 0.45 
)18j . In our model, fin = 87rG'p(M)/3'Uoj where uo is the present value of the Hubble parameter. 
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Consequently, the upper bound to \vo\ is (10~^^-10^^°) x yr~^, and the upper 
bound to \{v/u)o\ is about (0.01-1) x n~^h^^ (where we have used Uq = 100 h km 
s"^ Mpc~^). (For h = 0.71 as suggested by WMAP [21 and the special case with 
n = 3, the upper bound ranges from 0.0047 to 0.47.) Therefore the present state of 
the universe should be around the ka/ a^u^-sxis in the parameter space. 

The above constraint on \{v/u)o\, the requirement do > from SN la data, and 
0.15 ^ Qm ^ 0.45 [IH] together imply that the present state of the universe should be 
in the up-left direction of the fixed point at (—1,0) in the parameter space (Fig. IH)), 
and its evolution should follow type-(a) or type-(d) trajectories in Fig. HI Type-(d) 
trajectories are much more favored than type (a) for two reasons: Firstly, type-(d) 
trajectories possess smaller \v/u\ through this era so that they are more consistent 
with the bound to the variation of the Newtonian gravitational constant. Secondly, 
it has been shown in the previous discussions that in the radiation-dominated era 
the universe may approach to the states with stable extra dimensions (i.e. v/u ^ 
0), which should be the most probable initial conditions for the succeeding matter- 
dominated era, accordingly leading to type-(d) trajectories. 




Figure 5: An enlarged plot for the region where (d)-type trajectories pass. The 
trajectory (white dashed curve) which entails the largest acceleration parameter % is 
plotted. 

There is a boundary between type-(a) and type-(d) trajectories. The intersection 
of this boundary and the present condition 0.15 ^ ^{m) ^ 0.45 in the parameter space 
provides an upper bound, which is about 6.09 x 10~^, to the acceleration parameter 
go- These facts are summarized in Fig. |31 where the point (—1.173,0.165) is the 
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intersection of the curve fl(M) = 0.35 and the boundary between type-(a) and type- 
(d) trajectories, and the white dashed hne denotes one of type-(d) trajectories which 
pass through the neighborhood of this intersection point. In particular, we note 
that the maximal-acceleration-parameter state represented by this intersecting point 
entails small enough v/u {v/u = 0.165) so that a moderate bound to the variation of 
the Newtonian gravitational constant is satisfied. This fact indicates that the essential 
difficulty from the variation of the Newtonian gravitational constant can be somewhat 
avoided for type-(d) trajectories. Nevertheless, this value of Qq is smaller than that 
required by SN la data together with the results of CMB and LSS observations within 
the framework of the standard cosmology According to the result in jTHj, the 
most probable values of {Qm,u!x) corresponding to the requirement go = 6.09 x 10"^ 
is about (0.35, —0.51), which is roughly on the margin of the 2.5a confidence interval. 
(In [in] wx denotes the equation of state of dark energy.) 

Scale Factor 




■0.01 -0.00 



0.006 -0.004 -0.002 



Time 



Figure 6: The evolution of the scale factor a{t) for trajectories in the ED and the SC 
model. 

Now we want to compare the trajectory in Fig. which entails the maximal value 
of go among type-(d) trajectories in our model with extra dimensions, to be called 
"ED model" , with that predicted in the standard cosmology in which a fiat universe is 
considered to consist of pressureless matter with the density parameter Qm and dark 
energy with the density parameter fix (= 1 ~ ^m) and a constant state parameter 
Wx {wx = Px/px), to be called "SC model". In the ED model we need to specify the 
present condition about the essential parameters fl(^M), ^o, Y^o, and uq for obtaining 
an evolution path, i.e., a numerical solution of the Einstein equations (121)^© j while 
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in the SC model the values of Qm, wx, and Hq {Hq = uq by definition) are necessary 
information. In order to make these two evolution paths to be as similar as possible 
for the present time, we choose the present conditions in the following consistent 
way: First of all the acceleration parameter go in both models is specified to be 
6.09 X 10~^, i.e., the maximal value of go for type-(d) trajectories in the ED model. 
As mentioned in the last paragraph about the ED model, the intersection of the 
curve fl[M) = 0.35 and the (white dashed) trajectory in Fig. Elis at (—1.173,0.165), 
which provides the required consistent value of {Xq,Yq). In summary, we choose the 
following present conditions: fl(M) = 0.35, Xq = —1.173 and Yq = 0.165 for the ED 
model and Qm = 0.35 and Qx = —0.51 for the SC model, as well as Hq = uq = 71 
km s~^ Mpc~^ from WMAP Numerical solutions are shown in Fig. IHl These two 
evolution paths have similar behavior at the present time but diverge in the earlier 
time. In particular, the ED path entails a little shorter age. Possible observational 
signatures from these discrepancies in the age and the earlier-time behavior are under 
investigation. 

4 Conclusions and Discussions 

We have discussed the evolution of a universe under the infiuence of fiat extra dimen- 
sions in various eras via a well-chosen parameter space that provides a systematic, 
useful, and convenient way for analysis. Since we assume the fiatness (zero curvature) 
of the extra space, the only ingredient of extra dimensions that can affect the evolu- 
tion of the universe is the evolution of their size. We have studied general features of 
the evolution, thereby exploring whether the picture is viable, subject to constraints 
from observations, for describing our world. 

Through the above studies we have found a natural evolution pattern from the 
blazing era to the matter-dominated era in this model. In the blazing era, the ordinary 
space and the extra space tend to synchronize their expansion rates, i.e., v/u 1, 
and meanwhile the universe is decelerating, presuming that an extremely tiny curva- 
ture contribution ka/a^v^ has been guaranteed initially by infiation.^ This tendency 
provides a natural initial condition, {ka/a^u^ , v/u) = (±e, 1), for the succeeding 
radiation-dominated era (where '±e' indicates a tiny curvature contribution). 

In the radiation-dominated era, the expansion rate of the extra space, v, tends to 
approach zero from ti/w ~ 1 (as suggested above), meanwhile maintaining the decel- 
erating phase. ^ It is a good feature for building models with extra dimensions that 
the expansion rate of the extra space decreases to zero automatically, that is, extra 
dimensions are stabilized automatically in this era without resorting to any artificial 

•^This tendency corresponds to the fixed point at {ka/a?v? , v/u) ~ (0, 1) in the parameter space 
for tlie blazing era, which is stable in the Y direction but unstable in the X direction. 

''This tendency corresponds to the fixed point at (ka/a^u^ , v/u) = (0, 0) in the parameter space 
for the radiation-dominated era, which is stable in the Y direction but unstable in the X direction. 
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mechanism. In particular, for small enough v we can recover the standard cosmol- 
ogy (without extra dimensions) and the essential predictions therein (in particular, 
Big Bang Nucleosynthesis) in the radiation-dominated era. Therefore in our model 
with extra dimensions the predictions in the standard cosmology in the radiation- 
dominated era can be preserved in a natural way. 

In the matter-dominated era, a natural initial condition, [ka/a^u^ , v/u) = (±e, 0), 
is provided by the preceding, radiation-dominated era. For the case of a negative 
curvature, i.e. with the initial condition (— e, 0), the universe will eventually change 
its initial decelerating phase to the accelerating one and approach the attractor at 
{ka/a?v? , v/u) = (—1, 0) with stable extra dimension and significant negative curva- 
ture contribution. As we have shown, this evolution pattern is marginally consistent 
with the limit on the variation of the Newtonian gravitational constant. The maximal 
acceleration associated with this evolution pattern deviates from the one required by 
SN la data together with CMB and LLS observations (within the framework of the 
standard cosmology) [12] with 2.5a deviation. Thus this evolution pattern is still not 
ruled out by the present observations. 

We note that in this natural evolution pattern the expansion rate of extra dimen- 
sions will eventually approach zero, consequently extra dimensions being stabilized 
automatically without resorting to any artificial mechanism, in both the radiation- 
dominated and the matter-dominated era. This is the key ingredient for generating 
accelerating expansion of the present universe as well as satisfying the limit on the 
variation of the Newtonian gravitational constant (both caused by the evolution of 
extra dimensions), meanwhile preserving essential predictions in the standard cosmol- 
ogy for the radiation-dominated era. In particular, we emphasize that the naturalness 
of this evolution pattern, whose behavior is not sensitive to the initial conditions, in- 
dicates a solution to the cosmic coincidence problem ("why now" problem) of dark 
energy — why dark energy starts dominating the universe now, or, more precisely, 
why the accelerating phase starts at the present epoch. The existence of this natural 
evolution pattern implies that the late-time accelerating expansion of an open uni- 
verse ^ in the matter-dominated era is guaranteed, accordingly the cosmic coincidence 
problem being solved, in our flat ED model with the help of inflation in the very early 
time. 

Nevertheless, type-(d) trajectories in general entail a too short age of the universe 
and negative curvature in the ordinary space rather than almost zero curvature sug- 
gested by CMB data within the framework of standard cosmology. These features 
might eventually rule out this evolution pattern for describing our universe. Further 
detailed studies about the consistency with observations are in progress. Moreover, 
in order to complete the picture of the scenario for a universe with extra dimensions 
and find the best-fit trajectory therein, a more general model with nonzero curvature 
in the extra space is also under investigation. 

^Actually this open universe is nearly flat in the earlier time. 
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